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Polyakov loop correlations at finite temperature in two-flavor QCD are studied in lattice simu- 
lations with the RG-improved gluon action and the clover-improved Wilson quark action. From 
the simulations on a 16 3 x 4 lattice , we extract the free energies, the effective running coupling 
g e s(T) and the Debye screening mass mo(T) for various color channels of heavy quark— quark 
and quark— anti-quark pairs above the critical temperature. The free energies are well approxi- 
mated by the screened Coulomb form with the appropriate Casimir factors. The magnitude and 
the temperature dependence of the Debye mass are compared to those of the next-to-leading order 
thermal perturbation theory and to a phenomenological formula given in terms of g e ff(T). Also 
we made a comparison between our results with the Wilson quark and those with the staggered 
quark previously reported. 



XXIVth International Symposium on Lattice Field Theory 
July 23-28, 2006 
Tucson, Arizona, USA 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Static quark free energies at finite temperature with two flavors of improved Wilson quarks Y. Maezawa 



1. Introduction 

The interaction between heavy quarks in a hot QCD medium is one of the most important 
quantities to characterize the properties of the quark-gluon plasma (QGP). Experimentally, it is 
related to the fate of the charmoniums and bottomoniums in the QGP created in relativistic heavy 
ion collisions. In this report, we present our recent studies on the heavy-quark free energy in 
dynamical simulations of two-flavor QCD with the Wilson fermion. We study the free energy of 
a quark (Q) and an antiquark (Q) separated by the spatial distance r in the color singlet and octet 
channels, and also study the free energy of Q and Q in the color anti-triplet and sextet channels. We 
adopt the Coulomb gauge fixing for the gauge-non-singlet free energies. By fitting the numerical 
results with the screened Coulomb form, we extract an effective running coupling and the Debye 
screening mass in each channel as a function of temperature. 

We find that (i) the free energies in the different channels at high temperature (T <; 2T C ) can 
be well described by the channel-dependent Casimir factor together with the channel-independent 
running coupling g c s{T) and the Debye mass mr}(T), (ii) the next-to-leading order result of the 
Debye mass in thermal perturbation gives better agreement with mo(T) on the lattice than that 
of the leading order, (iii) g e a(T) and mjj(T) may be related through the leading order formula, 
mo{T) = Wl+Nf/6 ges(T)T, so that most of the higher order and non-perturbative effects on the 
Debye mass may be absorbed in g e s(T) for T > l.5T c , and (iv) there is a quantitative discrepancy 
of mj)(T) between our results using the Wilson quark and those using the staggered quark even at 
T ~ 4T C . 

2. Lattice action and simulation parameters 

We employ the renormalization group improved gluon action and clover improved Wilson 
quark action with two-flavors. The line of constant physics, on which the quark mass (or the ratio 
of pseudo-scalar and vector meson masses) is kept fixed in the space of the coupling j8 and the 
hopping parameter K, has been studied in refs. [jl|] using the same action and is discussed further 
in ref. We perform the simulations on the lines of constant physics for the quark masses 
corresponding to m n /m p = 0.65 and 0.80. Ten (seven) different temperatures are taken in the 
interval T = 1.00r c ~ 4.02T C for m % /m p = 0.65 (T = l.01T c ~ 3.01T C for m^/m p = 0.80) on a 
N s x N t = 16 3 x 4 lattice. The hybrid Monte Carlo algorithm is employed to generate full QCD 
configurations, and the free energy in each color channel is measured using 500 configurations 
at every ten trajectories after thermalization of 400-1000 trajectories. The statistical errors are 
determined by a jackknife method with the bin-size of 10 configurations. 

3. Heavy quark free energies 

The free energy of static quarks on a lattice may be described by the correlations of the 
Polyakov loop: £2(x) = n^Li^4( T ) x ) where N t is a lattice size in the temporal direction, and 
the Ufi(x,x) £ SU (3) is the link variable. By appropriate gauge fixing (such as the Coulomb gauge 
fixing), one can define the free energy in various color channels separately [fj]: the color singlet QQ 
channel (1), the color octet QQ channel (8), the color anti-triplet QQ channel (3*), and the color 
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Figure 1: Simulation results of the normalized free energies scaled by T c for color singlet and octet QQ 
channels (left) and color anti-triplet and sextet QQ channels (right) at m n /m p = 0.65 and several tempera- 
tures. 



sextet QQ channel (6) : 

e- F ^l T = i<TrQt(x)Q(y)), (3.1) 
e -Fk(r,T)/T = I( T rnt(x)Trn(y)) - l(Tr^(x)n(y)), (3.2) 
e -F 6 (r.T),T = _L (Tra(x)TrQ(y)) + l(Trn(x)n(y)), (3.3) 

e -F 3 ,(r,T)/T = L Tra ( x ) Tra (y)) _ i(Tra( x )n( y )), (3.4) 

6 6 

where r = |x — y|. 

We introduce normalized free energies (Vi,V$,V6,V3*) which are expected to approach zero at 
large distances above T c . This is equivalent to defining the free energies by dividing the right-hand 
sides of Eq. Q - (0) by (TrH) 2 . 



The normalized free energies are shown in Fig. [I] for color singlet and octet QQ channels 
(left) and color anti-triplet and sextet QQ channels (right) for m n /m p = 0.65 and T > T c . They 
are "attractive" in the color singlet and anti-triplet channels and "repulsive" in the color octet and 
sextet channels. Also they become weak at long distances as T increases due to the effect of Debye 
screening. These behaviors are qualitatively similar to the case in the quenched simulations with 
the Lorenz gauge reported in ref. [Q|. 

To study the screening effects in each color channel more closely, we fit the free energies by 
the screened Coulomb form: 

V M ^T)=C{M)^^e- m ^\ (3.5) 

where a e s(T) and mj)(T) are the effective running coupling and Debye screening mass, respec- 
tively. The C(M) = (£q = i t\ ■ t^M is the Casimir factor for each color channel (M) defined as 

C(l) = -i C{T) = - 2 -, C(8) = i C(6) = i. (3.6) 
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Figure 2: The effective running coupling Cfeff(T') (left) and Debye screening mass mo{T) (right) for each 
color channel as a function of temperature from the large distance behavior of the potentials at m K /m p = 
0.65. 



Here, it is worth stressing that, with the improved actions we adopt, the rotational symmetry is 
well restored in the heavy quark free energies [Q]. Therefore we do not need to introduce terms 
correcting lattice artifacts at short distances in eq. ( |3~3| ) to fit the data shown in Fig.l. 

The Debye screening effect is defined through the long distance behavior of Vm(i",T). In 
order to determine the appropriate fit range, we estimate the effective Debye mass from a ratio of 
normalized free energies: 



Mr;r) = -log^- T ^--log 



Ar 
1 + — 
r 



(3.7) 



Investigating the plateau of mj)(T\ r), we choose the fit range to be v 1 1/4 < rT < 1 .5. Systematic 
errors due to the difference of the fit range are about 10% for T ^ 2T C . 

The results of the a e ff(T) an d m D {T) are shown in Fig. ^|for m % jm p = 0.65. Similar behavior 
in both results are obtained for m n /m p = 0.80. We find that there is no significant channel depen- 
dence of OCes(T) and m D (T) at sufficiently high temperature (T J> 2T C ). In other words, the channel 
dependence of the free energy at high temperature may be well absorbed in the kinematical Casimir 
factor as first indicated in a quenched study M]. 



4. Debye mass on the lattice and that in perturbative theory 

Let us first compare the Debye mass on the lattice with that calculated in thermal perturbation 
theory. First of all, the 2-loop running coupling is given by 

g 2l V) = ftln 2 + | Into ( JL) 2 , ( 4 . 1) 

where the argument in the logarithms may be written as fi/Aj^ = (jx/T)(T /T C )(T C /Aj^) where 
we adopt T c /A^ 2 ~ 0.656 for the last factor |[T[ |6|]. The renormalization point /I is assume to be 
in a range /I = %T — 2>nT . By using g2\ as a function of T /T c , the Debye screening mass in the 
leading-order (LO) thermal perturbation is given as m^°(r)/r = Nf/ 6 gn(T), where the 
effect of the quark mass is neglected. 
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Figure 3: The Debye screening masses mo(T) at m K /m p = 0.65 and 0.80 in the color singlet channel 
together with that calculated in the leading-order (left) and next-to-leading-order (right) thermal perturbation 
theory shown by the dashed lines, pi is the renormalization point chosen at ji = nT, 2kT, 3nT. 



Figure ||](left) shows the moiT) in the color singlet channel compared with m^°(r)/r for 
/I = nT, 2nT and 3nT. We find that the screening mass m]f(T) in the leading order perturbation 
theory does not reproduce the lattice data, which has been known in the quenched QCD [f7|] and in 
the full QCD with the staggered quark []8|] . 

To study higher-order contributions in the thermal perturbation theory, we consider the Debye 
mass in the next-to-leading-order calculated by the hard thermal resummation given in ref. [0J. 



m 
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111 



mag 



+ 



(4.2) 



Here m mag denotes the magnetic screening mass assumed to be of the form m mag (T) = C m g 2 (T)T. 
Since the factor C m cannot be determined in the perturbation theory due to the infrared problem, 
we adopt C m ~ 0.482 calculated in quenched lattice simulations [Q] as a typical value. ( If we fit C m 
from our lattice data m D (T = 4.02r c .) with ju = 2kT, we obtain C m ~ 0.40). m^ LO (T) for different 
choice of the renormalization point is shown in Fig. ||(right) by the dashed lines together with the 
lattice data. They have approximately 50 % enhancement from the leading order results and lead 
to a better agreement with the lattice data. 



5. Phenomenological relation between a e ff and mo 

So far, we have fitted the free energies on the lattice with a e ff and m D as independent parame- 
ters. Here let us introduce an "effective" running coupling as g e s(T) = yAnot^T). Suppose that 
mj){T) is expressed by g e ff(T) according to the leading-order perturbation, 



(5.1) 



This relation means that the following ratio R evaluated from our lattice data should be close to 
unity, i.e. R(T) = (1 +A^/6)~ 1 / 2 {m D {T) /T) / ^ 47ia sS (T) ~ 1. 
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Figure 4: The ratio R(T) in the text which is supposed to be close to unity if mo{T) = y/l + N//6 g e s(T)T 
holds. The data are for the color singlet channel. 

In Fig. ||, the lattice data of R(T) for the singlet channel are shown as a function of T. We 
find that R(T) is consistent with unity even at T <; l.5T c with 10% accuracy. This is a non-trivial 
observation particularly near T c and suggests that the major part of the higher-order effects and 
non-perturbative effects of niD{T) can be expressed by the effective running coupling g e s(T). We 
note that a similar effective coupling defined through the lattice potential was discussed to improve 
the lattice perturbation theory at T = [|T0|]. 

6. Comparison with the staggered quark 

Finally, we compare the results of a e ff (T) and m D (T) obtained with the Wilson quark action 
(present work) with those with an improved staggered quark action. The latter simulation was 
done on a 16 3 x 4 and with a quark mass corresponding to m n /m p ~ 0.70 [§]. The comparison 
is shown in Fig. || for a e ff(T) (left panel) and mjj{T) (right panel). Although a e s(T) does not 
show significant difference between the two actions, mo(T) in the Wilson action is systematically 
higher than that of the staggered action by about 20% even at 4T C . This difference should be further 
investigated by increasing the temporal lattice size. 

7. Summary 

We have studied the free energy of QQ and QQ systems in 2-flavor QCD at finite temperature 
using the lattice simulation with the renormalization group improved gluon action and the clover 
improved Wilson quark action on a 16 x 4 lattice. The free energy normalized to be zero at large 
separation show attraction (repulsion) in the color singlet and anti-triplet channels (color octet and 
sextet channels). 

The screened Coulomb form with the Casimir factor and with the effective coupling a e ff(T) 
and the Debye screening mass mr>{T) as free parameters is used to fit the free energy in each 
channel. (Xeff(T) and moiT) become universal and all the channel dependence is absorbed in the 
Casimir factor for T <; 2T C . The magnitude and the T-dependence of the Debye mass moiT) 
is consistent with the next-to-leading order calculation in the perturbation theory as well as the 
leading order perturbation with the "effective" running coupling. 
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Figure 5: Comparison of the 0feff(T)(left) and mo(7')(right) between the results of the Wilson quark action 
and staggered quark action. 



The results from an improved Wilson quark action are compared with these from an improved 
staggered quark action with the same lattice size and similar quark mass. The a e f{(T) does not show 
appreciable difference between the two actions, whereas the wd(T) of the Wilson quark action is 
systematically higher than that of the staggered quark action by 20%. The simulations with larger 
lattice sizes especially in the temporal direction (such as N t = 6 and larger) should be carried out 
as well as those at small quark masses. 

Acknowledgements: We thank O. Kaczmarek for providing us the data from a staggered quark 
action. This work is in part supported by Grants-in-Aid of the Japanese Ministry of Educa- 
tion, Culture, Sports, Science and Technology, (Nos. 13135204, 15540251, 17340066, 18540253, 
18740134). SE is supported by the Sumitomo Foundation (No. 050408), and YM is supported by 
JSPS. This work is in part supported also by the Large-Scale Numerical Simulation Projects of 
ACCC, Univ. of Tsukuba. 



References 

[1] A. Ali Khan et al. (CP-PACS Collaboration), Phys. Rev. D63 (2001) 034502; 64 (2001) 074510. 

[2] N. Ukita, S. Ejiri, T. Hatsuda, N. Ishii, Y. Maezawa, S. Aoki and K. Kanaya, in proceedings of Lattice 



2006. 



PoS (LAT2006) 150 



[3] S. Nadkarni, Phys. Rev. D33 (1986) 3738; 34 (1986) 3904. 

[4] A. Nakamura and T. Saito, Prog. Theor. Phys.lU (2004) 733; 112 (2004) 183. 

[5] S. Aoki et al. (CP-PACS Collaboration), Phys. Rev. D60 (1999) 1 14508. 

[6] M. Gockeler et al, Phys. Rev. D73 (2006) 014513. 

[7] A. Nakamura, T. Saito and S. Sakai, Phys. Rev. D69 (2004) 014506. 

[8] O. Kaczmarek and F. Zantow, Phys. Rev. D71 (2005) 1 14510. 

[9] A. K. Rebhan, Phys. Rev. D48 (1993) 3967. 

[10] G. P. Lepage and P. B. Mackenzie, Phys. Rev. D4S (1993) 2250. 



7 



